A 62 Day X-Ray Periodicity and an X-Ray Flare from the 
Ultraluminous X-Ray Source in M82 
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ABSTRACT 

In 240 days of X-ray monitoring of M82, we have discovered an X-ray periodicity at 62.0 ±2.5 
days with a peak to peak amphtude corresponding to an isotropic luminosity of 2.4 x 10**^ ergs"^ 
in M82 and an X-ray flare reaching a peak luminosity of 9.8 x 10^" ergs~^. The periodicity and 
flare likely originate from the ultraluminous X-ray source (ULX) in M82 which has been identifled 
as a possible intermediate mass black hole. We suggest that the 62 day modulation is due to 
orbital motion within an X-ray binary with a Roche-lobe overflowing companion star which would 
imply that the average density of the companion star is near 5 x lO^^gcm"'^ and is therefore a 
giant or supergiant. Chandra observations just after the flare show an energy spectrum that is 
consistent with a power-law with no evidence of a thermal component or line emission. Radio 
observations made with the VLA during the flare allow us to rule out a blazar identification for the 
source and place strong constraints on rclativistically beamed models of the X-ray emission. The 
Chandra observations reveal a second X-ray source reached a flux of 4.4 x 10~^^ ergcm"^ s~^ in 
the 0.3-7 keV band which is dramatically higher than any flux previously seen from this source and 
corresponds to an isotropic luminosity of 1.1 x 10"*" erg s~^. This source is a second ultraluminous 
X-ray source in M82 and may give rise to the QPOs detected from the central region of M82. 



Subject headings: black hole physics 
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The bright X-ray sources in external galaxies, 
known as ultraluminous X-ray sources (ULXs), 
are of considerable current interest because they 
may be 'intermediate-mass' black holes (Colbert 
& Mushotzky 1999; Makishima et al. 2000; Kaaret 
et al. 2001). One of the most extreme ULXs 
is located in the nearby starburst galaxy M82 
(Ptak & Griffiths 1999). Chandra observations 
showed that the brightest source in the galaxy, 
CXOU J095550.2-h694047 = X41.4-h60 (Kaaret et 
al. 2001; Matsumoto et al. 2001), is offset from 
the nucleus, not coincident with any radio super- 
nova remnant, and highly variable on time scales 
of weeks. The offset of the source from the dynam- 
ical center of the galaxy places an upper bound on 
its mass of 10^ — 10^ Mq, depending on its age, due 
to dynamical friction (Kaaret et al. 2001). These 
observations rule out a low luminosity AGN or a 



very luminous supernova or hypernova remnant as 
the source of the X-ray emission. If the radiation 
from X41. 4-1-60 is isotropic and from an accreting 
object, then the mass of the accretor must be in 
excess of 5OOM0. Combined with the mass upper 
bound, this would place the object firmly in the 
"intermediate-mass" black hole category. 

However, alternative explanations must be con- 
sidered, particularly that the ULXs radiate at 
super-Eddington luminosities Begelman (2002) or 
that the radiation is beamed, so the true luminos- 
ity is significantly less than calculated assuming 
isotropic emission. King et al. (2001) have sug- 
gested that ULXs are high-mass X-ray binaries in 
a short-lived phase which occurs early in the life 
of almost every system and during which mass is 
transferred from the companion to the compact 
object on the thermal time scale of the compan- 
ion - leading to super-Eddington accretion rates. 



The X-ray emission in this phase is mechanically 
beamed, producing high observed X-ray fluxes for 
observers near the beaming axis. However, the 
flux of X41. 4-1-60 is so high that producing it via 
mechanical beaming of a stellar-mass black hole is 
problematic. King et al. (2001) specifically point 
to X41. 4-1-60 as an example of a system where their 
model may not apply. 

Mirabel & Rodriguez (1999) pointed out that 
radio microquasars with their axes aimed at us 
( "microblazars" ) should occur in nearby galax- 
ies. Kording et al. (2002) have suggested that the 
ULXs may be stellar-mass black holes in which a 
relativistic jet produces most of the observed X- 
ray flux. Jets aligned nearly along our line of sight 
could then produce high apparent X-ray fluxes 
via relativistic beaming. If the ULXs are microb- 
lazars, then they should be radio sources and the 
X-ray and radio emission of the ULXs should be 
well correlated and highly variable on time scales 
as short as days. Kaaret et al. (2001) noted that 
an unusual radio transient, 41.5+59.7 (Kronbcrg 
& Wilkinson 1985), lies within the error box for 
X41. 4-1-60. This transient was detected in Febru- 
ary 1981 at a 5 GHz flux of 7 mJy and not de- 
tected with an upper limit of 1.5 mJy in October 
1983. The sharp decrease in flux indicates that 
it is not a supernova remnant. Transient radio 
emission from an accreting compact object would 
suggest the presence of a jet of relativistic particles 
(Kaaret et al. 2003). 

Recently, Strohmayer & Mushotzky (2003) 
found quasipcriodic oscillations (QPOs) in the 
range 50-110 niHz which they identify as arising 
from X41.4+60. If the QPO frequency is limited 
by the orbital frequency around a non-rotating 
black hole, then it would imply a mass limit of 
2 X lO"' Mq. The QPOs were discovered in XMM- 
Newton data and confirmed in RXTE data. The 
QPOs are only occasionally detected in the RXTE 
data and there is no apparent correlation between 
QPO detections and the source flux level. Based 
on a comparison with the spectral and timing 
properties of stellar-mass black hole X-ray bina- 
ries, Fiorito & Titarchuk (2004) estimate a mass of 
the order of lOOOM© for the compact object pro- 
ducing the QPOs. The QPOs were confirmed in a 
longer XMM-Ncwton observation, but found at a 
somewhat higher frequency of 113 mHz (Muccia- 
relli et al. 2006; Dewangan, Titarchuk, & Griffiths 



2006). 

The position of X41.4-(-60 is within 1" of the po- 
sition of the infrared source and super star cluster 
MGG 11 (Kaaret et al. 2004; Portegics Zwart et 
al. 2004a) . Simulations of the dynamical evolution 
of the cluster MGG 11 show that stellar collisions 
in its extremely dense core may have led to nu- 
merous stellar collisions and the formation of an 
intermediate mass black hole (Portcgies Zwart et 
al. 2004a). 

In order to obtain more information about 
X41.4+60, we monitored M82 for 240 days us- 
ing the Proportional Counter Array (PGA) on the 
Rossi X-Ray Timing Explorer (RXTE) with obser- 
vations obtained typically every 2 days. The PGA 
is not an imaging instrument, but is adequate to 
measure the X-ray flux and rapid variability of the 
integrated output of the entire galaxy. We discov- 
ered an X-ray periodicity from M82 at 62.0 ± 2.5 
days which we have previously reported (Kaaret, 
Simet, & Lang 2006). Additional details of the 
analysis are presented here. We also found a tran- 
sient event where the X-ray flux of the galaxy 
increased by a factor of three over the quiescent 
level over a period of a few days. After the onset 
of this remarkable X-ray flare, we obtained imag- 
ing X-ray observations with the Ghandra X-Ray 
Observatory and radio observations with the Very 
Large Array (VLA). We describe X-ray observa- 
tions with RXTE and Ghandra in § 2, the VLA 
data in § 3, and discuss the results in § 4. 

2. RXTE observations 

We obtained a total of 141 observations of M82 
with a total effective exposure of 281.2 ks using the 
Proportional Gounter Array (PGA) on the Rossi 
X-Ray Timing Explorer (RXTE) under program 
90121 (PI Kaaret). Two of the observations were 
made early to test our procedures for the acquisi- 
tion and analysis of realtime data. The bulk of the 
observations consisted of monitoring the galaxy 
typically every second day during the period over 
which the VLA was in the A, BnA, or B config- 
uration during 2004 and 2005 and covered dates 
from MJD 53252 to 53490. 

To enable target of opportunity observations 
with Ghandra and the VLA, we analyzed the real 
time data from RXTE soon after the data be- 
came available on the RXTE ftp site. The anal- 



ysis was carried out with the FTOOLS package 
(version 5.3.1) which is part of the HEAsoft soft- 
ware. After the data were retrieved, they were 
filtered to select good time intervals such that 
Proportional Counter Unit (PCU) 2 was on, the 
source was more than ten degrees above the hori- 
zon, the pointing offset from the source was less 
than 0.01°, the satellite was out of the South At- 
lantic Anomaly (SAA), and the electron contam- 
ination was less than 0.1. A background file was 
made for the same interval as the observation with 
the model pca_bkgd_cmfaintl7_eMv20031123.mdl 
and the newest predicted SAA passage file down- 
loaded at the same time as the data. The FTOOL 
saextrct was used to make both light curves and 
spectra for the observation and the background, 
using only data from PCU 2. The spectra were 
fitted with a power-law model using the spectral 
fitting program XSPEC, which is part of HEA- 
soft, to obtain a measurement of the flux in the 
2-10 keV band. Our flux monitoring led to the 
detection of an X-ray flare reaching a flux of 
4.5 X 10~^^ ergcm"^ s~^ in the 2-10 keV band on 
26 Jan 2005 (MJD 53396.4). On the basis of this 
flare detection, we triggered target of opportunity 
observations with Chandra and the VLA. 

After the observing program was completed, 
we re-examined all of the observations using the 
RXTE production data. This had the advantages 
that some data gaps in the real-time data were 
fllled and, more importantly, that the SAA pas- 
sage history, rather than the SAA passage predic- 
tions, could be used in estimating the PCA back- 
ground. To produce an energy spectrum for each 
observation, we used the Rex script ^. We pro- 
duced spectra using only the top layer in PCU 
2, including channels 2-44 which cover the energy 
range 2-20 keV. Response matrices were made for 
each observation. The expected background was 
estimated using the same model mentioned above 
and a systematic uncertainty was included in the 
background estimate ^. The spectrum was fitted 
using XSPEC with a power-law model with inter- 
stellar absorption with the column density fixed 
to 3 X 10^^ cm~^. The flux was calculated for the 
2-10 keV range. 
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Fig. 1.— Light curves of M82 obtained the PCA 
on RXTE. The plot shows the flux in the 2-10 keV 
band calculated for each observation versus the ob- 
servation date in MJD. The diamonds indicate the 
observations where QPOs were detected. The ver- 
tical dashed line indicates the time of the Chandra 
observation. The solid horizontal line indicates 
the interval over which VLA observations were ob- 
tained. A strong modulation in the X-ray flux near 
a period of 62 days is apparent. An X-ray flare oc- 
curred around MJD 53400. 

Fig. 1 shows the flux in the 2-10 keV band ver- 
sus time. The X-ray flare is obvious as the high 
flux points around MJD 53400. The peak flux 
occurred on MJD 53400.2 and reached a level of 
6.3 X lO^^^ergcm^^ s^^ in the 2-10 keV band. 
The duration of the flare is rather short, only a 
few days. This is in contrast to the light curve 
of M82 presented by Rephaeh & Gruber (2002), 
where the duration of the high flux interval is at 
least 100 days. The photon index is generally be- 
tween 2.2 and 2.5 and shows no strong correlation 
with flux. 

The light curve shows an apparent modulation 
with an amplitude of 1.6 x 10^^^ ergcm^^ s^^ and 
a period near 60 days. We calculated a peri- 
odogram, shown in Fig. 2, according to the method 
of Home & Baliunas (1986) with the power nor- 
malized by the total variance of the data. For 
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Fig. 2. — Periodogram calculated from the flux 
measurements of M82 obtained for individual ob- 
servations. A strong peak is apparent at a period 
of 62.0 ± 2.5 days. The periodogram is calculated 
via the method of Home & Baliunas (1986) with 
the power normalized by the total variance of the 
data. 

the periodogram, we used our entire RXTE data 
set. There is a peak at a period 62.0 days with 
a power of 27.9. The probability of chance oc- 
currence calculated assuming a Poisson distribu- 
tion, or white noise background, and taking into 
account the number of trials is 1.1 x 10^^°, cor- 
responding to a 6.7(7 detection. We estimate the 
90% confidence error on the period to be 2.5 days. 
There is a secondary peak at 125.7 days with a 
power of 15.3. Removing the two early observa- 
tions does not significantly affect the periodogram. 
The peak is then at 61.6 ±2.7 days with a power of 
27.5. Removing the observations during the X-ray 
flare near MJD 53400 does not shift the peak, but 
increases its power. With the points with fluxes 
above 4 x 10~^^ erg cm~^ s~^ removed, the peak is 
at 62.0 ± 2.1 days with a power of 32.8. Removal 
of the two early points and the high flux points 
leads to a peak is at 61.6 ± 2.2 days with a power 
of 32.8. 

Accreting compact objects often exhibit red 
noise, or increasing timing noise at decreasing fre- 



quency with a power-law dependence of power on 
frequency. Red noise can produce apparent peri- 
odicities at periods comparable to the observation 
duration (Israel & Stella 1996). We tested the sig- 
nificance of the observed signal using a red noise 
background. First, we fit a power-law to the power 
versus frequency relation for periods shorter than 
40 days (Vaughan 2005). We found a spectral in- 
dex of -0.78. We then generated red noise with 
this spectral index and with a mean and variance 
equal to the mean and variance of the actual data. 
This procedure is conservative because it includes 
the variance due to the signal itself and also due 
to the X-ray flare. We generated simulated light 
curves using the rndpwrlc routine of the aitlib 
IDL subroutine library-^ provided by the Institut 
fiir Astronomie und Astrophysik of the Univer- 
sitat Tiibingen. For each light curve, we generated 
2048 data points with uniform spacing equal to one 
third of the average spacing of the main 139 obser- 
vations of our monitoring program and used only 
141 points from the middle of the simulated data 
with relative times matching the actual observa- 
tions. Generating a simulated light curve longer 
than the portion used helps eliminate the effects 
of red noise leakage. The simulated data were pro- 
cessed with the same procedures used to analyze 
the real data. We generated 100,000 trial light 
curves and searched for cases where the power at 
a period of 80 days or less, corresponding to three 
oscillation cycles in the main portion of the mon- 
itoring, was greater than or equal to the observed 
value of 27.9. We found only 2 cases, equivalent to 
a probability of chance occurrence of our observed 
signal of 2 X 10~^. 

The Q value of the peak, the period of the peak 
divided by the full width at half maximum power, 
in the periodogram is 4.5, which is fully consistent 
with that expected for a periodic process given 
the observation duration, but is only weakly con- 
straining. To attempt to determine whether or 
not the signal is strictly periodic, we measured 
the arrival times of the peak maxima. We flt the 
data near each peak, ±20 days, with a second or- 
der polynomial. We excluded the very high flux 
points with fluxes above 4 x 10^^^ ergcm^^ s~^. 
Comparing the peak times to an assumed 62 day 
period, we find that the measured arrival times 



^available at http://astro.uni-tuebingen.de/software/aitlib 



of first, second, and fourth peaks are consistent 
with the predicted arrival time for periodic sig- 
nal within 2 days. The third peak differs by 4 
days. However, the polynomial fit is affected by 
two high flux points in the decaying part of that 
peak (at MJD 53408.2 and 53412.7). If we fit the 
data within ±11 days of that peak, we find that 
the peak arrival time agrees within 1 day of the 
time predicted for a periodic signal. We note that 
we had observations only every second day and 
there are variations in the flux from other sources 
in M82 which are not related to the 62 day peri- 
odicity and which lead to uncertainty in the time 
of arrival determination. To estimate the uncer- 
tainty in the polynomial fits, we varied the width 
of the data windows and the placement of the data 
windows, i.e. offsets the windows with respect to 
the peak. We conclude that the peak arrival times 
are accurate to about 2 days. Therefore, the ar- 
rival times of the peaks are consistent with being 
periodic at an interval of 62.0 days. However, ad- 
ditional observations are essential to test the peri- 
odic nature of the signal. 

There is archival RXTE data from a previous 
set of observations of M82 made in 1997. How- 
ever, the time coverage of that data set is uneven 
and generally less dense than in our monitoring 
program. There is a hint of one peak in the 1997 
light curve around MJD 50585 which has a shape 
similar to the four seen in our data. However, the 
1997 data alone is insufficient to con&m or deny 
the periodicity. If we make a periodogram com- 
bining the two data sets, we find a peak at 61.5 
days. 

To examine the rapid variability of the source, 
152 event files were examined individually accord- 
ing to the following procedure using a program 
written in the Interactive Data Language (IDL). 
Events in the 2-10 keV energy band were selected 
and split into segments of 256 s. An FFT with 
a time resolution of 2~^ s was calculated for each 
segment. The FFTs were added incoherently. The 
resulting total power spectrum was logarithmi- 
cally rebinned with the bin width adjusted on a 
case by case basis. The IDL routine Imfit was 
then used to fit a model consisting of an expo- 
nential background and, if necessary, a Lorentzian 
peak to model a QPO. 

Of 152 event files, 10 contained QPOs with 
a single trial significance of Aa or greater. The 
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Fig. 3. — Selected power spectra of M82 obtained 
with the PCA on RXTE. Each power spectrum 
is labeled with the MJD of the observation date. 
The power spectra shown are those with QPOs 
detected with a single trial significance of Aa or 
greater. The QPO frequencies lie in the range 
from 80 to 190 mHz. 



power spectra from observations with QPO de- 
tections are plotted in Fig. 3. The centroid fre- 
quencies of the QPOs range from 80 to 190 mHz. 
The maximum modulated flux in the QPO is 
1.9 X 10~^^ergcm^^s^^ in the 2-10 keV band. 

2.1. Chandra observations 

After detection of the large X-ray flare in the 
PCA monitoring, a Chandra Target of Oppor- 
tunity Observation was triggered. The request 
was sent to the Chandra X-ray Center on 26 
Jan 2005 and the observation (ObsID 6097; PI 
Kaaret) began on 4 Feb 2005 at 23:34:40 UT (MJD 



53405.982). The observation was made using 
the Advanced CCD Imaging Spectrometer spec- 
troscopy array (ACIS-S; Bautz et al. 1998) and the 
High Resolution Mirror Assembly (HRMA). We 
obtained 52,773 s of useful exposure. The ACIS-S 
was used in imaging mode. The high flux antic- 
ipated from the source would have caused severe 
pile-up if the source were placed on-axis and the 
standard observing mode were used. To reduce the 
effect of pile-up, the target was offset 3'. 57 along 
the Y-detector coordinate. This placed the target 
off axis where the point spread function broadens 
and the X-ray flux is spread over multiple pixels in 
the CCD detector which reduces pile-up. In addi- 
tion, we operated only the S3 chip and employed 
a 1/8 sub-array mode to reduce the frame time to 
0.441 s and further reduce the pile-up. To check 
if pile-up was significant in the data acquired, we 
calculated the count rate in the 3x3 pixel cell 
at the peak of the brightest source in the image. 
The rate was 0.18 counts/frame, giving a pile-up 
fraction of 3.2% within this 3x3 pixel cell. The 
counts in the peak 3x3 pixel cell are 70% of the 
total source counts, so the overall pile- up fraction 
is 2.2%. Therefore, the effect of pile-up on the 
timing and power spectra is small. 

The Chandra data were subjected to standard 
data processing and event screening (ASCDS ver- 
sion 7.5.0.1 using CALDB version 3.0.0). The to- 
tal rate on the S3 chip was between 2.5-3.8 c/s 
for the entire observation and did not show any 
strong background flares. We constructed images 
for the 0.3-8 kcV and 2-8 keV bands using all 
valid events and used the wavdetect tool in the 
CL40 version 3.1 Chandrda data analysis package 
to search for X-ray sources. We found that use of 
the higher energy band seemed to produce more 
reliable point source detections and that the dif- 
fuse emission had less influence on the calculated 
source positions. Therefore, we used the source 
positions found from the 2-8 keV data in the re- 
mainder of the analysis. 

Fig. 4 shows the image in the 0.3-8 keV band. 
Sources with detection significance of 7a or higher 
are marked. The point sources appear as ellipses 
due to the off-axis pointing of the observation. 
The circle in the figure represents the extraction 
region used by Strohmayer & Mushotzky (2003) in 
their analysis of the XMM-Newton data. Several 
bright X-ray sources fall in or near the extraction 
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Fig. 5. — Light curves of the two brightest sources 
detected in the Chandra data. The upper curve 
is for the brightest source, X41.4+60. It shows 
a gradual increase on time scales of tens of kilo- 
seconds. The lower light curve is for the second 
brightest source, X42.3+59. It shows variability 
on time scales of kiloseconds. The light curves in- 
clude events in the 0.3-10 keV band. The time 
bins are 450 s. 



region and must have contributed to the energy 
and timing spectra attributed to X41. 4-1-60. 

The two brightest sources have fluxes which 
are an order of magnitude greater than any other 
source in the fleld. They are of particular in- 
terest and are the topic of the reminder of our 
Chandra analysis. The brightest source is CXOU 
J095550.2-h694047, the known ULX in M82 which 
was discussed in the introduction. Following 
Kaaret et al. (2001), we will refer to this source 
as X41.4+60. The second source is located at 
of a = 09h 55m 51s.040 and S ^ +69° 40' 
45".49 (J2000). This source is hkcly the same 
as CXOM82 095551.4-^694045 in Matsumoto et 
al. (2001) and may be the same as CXOM82 
095551.4+694043 in Grifliths et al. (2000). Fol- 
lowing the convention of naming sources in M82 
by their offset from a = 09h 55m 00s, S = +69° 40' 
00" (B1950), we refer to this source as X42.3+59. 
We note that Strohmayer & Mushotzky (2003) 




Fig. 4. — X-ray image of M82. The image includes all events in the 0.3-8 keV band. The light blue crosses 
indicate the positions of point sources. The brightest source is X41.4-I-60 and is located at the center of 
the image and labeled source # 1. The second brightest source is X42.3-I-59 is labeled source # 2. The 
light blue circle indicates the extraction region used for the analysis of XMM-Newton data by Strohmayer & 
Mushotzky (2003). Note that both of the brightest sources and several additional weaker sources lie within 
that extraction region. The arrow points north and has a length of 10". 



refer to the brightest source in M82 as M82 X- 
1. They assumed that this source is the same 
as X41.4-I-60 and ascribed it to be the origin 
of the QPOs found in the XMM-Newton and 
RXTE data. However, as discussed further be- 
low, X42. 3-1-59 lies within the extraction radius 
used by Strohmayer & Mushotzky (2003) to de- 
fine the X-ray counts from M82 X-1. Therefore, 
their source M82 X-1 may include contributions 
from both X41.4-h60 and X42. 3-1-59. 

We extracted photon lists and created light 
curves and energy spectra for both sources using 
standard CIAO tools. The light curves are shown 
in Fig. 5. The light curves represent the count rate 
in the 0.3-10 keV band. The upper curve is for the 
brighter source, X41. 4-1-60. The source appears to 



vary gradually on time scales of several thousand 
seconds. The flux of X42. 3-1-59 appears to have 
variations on faster times scales, but is relatively 
constant on long time scales. 

Fig 6 shows the timing power spectra for both 
sources derived from counts in the 2-10 keV band. 
The plotted spectrum for each source is the sum of 
16 individual power spectra with 8192 points each. 
The summed power spectrum was rebinned into 
logarithmically spaced bins before plotting. The 
sampling time was equal to the ACIS frame time 
of 0.441 s. In the power spectra plot, the curve for 
X41. 4-1-60 is the lower one. The timing noise from 
X41. 4-1-60 is consistent with Poisson fluctuations 
at frequencies from 0.001 Hz to 1 Hz. In contrast, 
X42. 3-1-59 shows significant noise near 1 niHz and 
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Fig. 6. — Power spectra of the two brightest 
sources detected in the Chandra data. The lower 
curve is for the brightest source, X41.4+60. Its 
power spectrum is consistent with Poisson noise 
over the frequency range plotted. The upper curve 
is for the second brightest source, X42.3+59. It 
shows significant timing noise at low frequencies 
and a possible QPO near 70 niHz. The power 
spectra are calculated from events in the 2-10 keV 
band. 

a possible QPO near 70 mHz. Taking into account 
the number of trials, i.e. the number of frequency 
bins in the power spectrum, the QPO at 70 mHz 
has a chance probability of occurrence of 9 x 10""^. 
The energy spectrum for X41.4+60 is shown in 
Fig. 7. We fitted the spectrum using the Sherpa 
fitting package which is part of CIAO. The spec- 
trum is adequately fitted by a power-law with 
interstellar absorption except for a region where 
there are significant residuals just above 2.1 keV. 
There is a strong feature in the instrument re- 
sponse at 2.1 keV corresponding to the Ir M edge. 
An analysis of grating data presented by Herman 
Marshall at the 2003 Chandra Calibration Work- 
shop^ shows that a r^ 10% correction is needed in 
the HRMA response at energies just above the Ir- 
M edge. Since this correction may depend on the 
off-axis angle of the source and does not appear 
to be adequately modeled in the current response 
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Fig. 7. — Energy spectrum of X41. 4-1-60. The 
spectrum is adequately fitted by a power-law with 
interstellar absorption except for a region just 
above the Ir M edge at 2.1 keV where the response 
matrix at large off-axis angles is not well deter- 
mined. 



matrices, we have chosen to remove the data for 
energies from 2.1-2.7 keV in performing the fits. 

For X41. 4-1-60, we find that the energy spec- 
trum is adequately fitted with a power-law with 
interstellar absorption giving x^/DoF = 72.8/66. 
The best fit parameters were a photon index 
of r = 1.67 ± 0.02 and an equivalent hy- 
drogen absorption column density of Nh ~ 
1.12 ± 0.02 X 10^2 cm-2 ^he source fiux was 
1.3 X lO^^ergcm-^s"^ in the 0.3-7 keV band 
and 1.5 x 10^^^ erg cm^^ s~^ in the 2-10 keV band. 
We note that an absorbed multicolor disk black- 
body model is strongly rejected with x^/DoF = 
4.0. Addition of a multicolor disk blackbody to 
the power-law model does not significantly im- 
prove the fit. For disk temperatures greater than 
0.5 keV, the normalization of the disk blackbody 
component is A'^ < 1. However, due to the signif- 
icant absorption at low energies, we cannot rule 
out the possible presence of a cool disk component 
with a temperature below 0.3 keV. Addition of an 
Fe-K emission line also does not improve the fit. 
Taking a gaussian line profile with a centroid of 



6.55 keV and a width of 0.33 keV to match the 
line reported by Strohmayer & Mushotzky (2003) , 
we place a 99% confidence level upper bound of 
70 eV on the equivalent width of the line. This 
is below the equivalent widths found from the 
XMM-Newton data. 

For X42.3+59, the source flux was 4.4 x 
10"^^ ergcm"^ s~^ in the 0.3-7 keV band and 
6.9 X 10~^^ ergcm"^ s~^ in the 2-10 keV band. A 
single power-law with interstellar absorption did 
not give an adequate fit. We found an adequate 
fit with the sum of two power-laws with indepen- 
dent interstellar absorption, with Fi = 1.33±0.07, 
{Nh)i = 2.8 ± 0.2 X 10^2 cm-^ F2 ^ 4.13 ± 0.01, 
Inh)2 = 9.0 ± 0.9 X 10^1 cm-2, and xV^oF = 
156.0/162. We also found an adequate fit with 
the absorbed sum of a multicolor disk blackbody 
plus a power-law. The best fit parameters were 
Nh = 2.29 ± 0.11 X 10^2 cm-^ F = 1.19 ± 0.06, 
disk inner temperature kTi„ = 0.072 ± 0.002 kcV, 
disk normahzation A^ = 3.2 x 10^, and x^/^oF = 
162.1/163. This disk normalization would imply 
a black hole mass in excess of 4 x 1O'''M0. 

3. VLA Observations 

Four separate radio observations using the Very 
Large Array (VLA) of the National Radio Astron- 
omy Observatory (NRAO)^ were made of M82 af- 
ter detection of X-ray fiarc. The first observation 
occurred on 29 Jan 2005 (MJD 53399) and the last 
on 5 Feb 2005 (MJD 53406). The radio observa- 
tions began one day before the peak X-ray flux and 
continued to 6 days after. All observations were 
made at 8.5 GHz, with integration times ranging 
from 1 to 2.5 hours and the array was in its hy- 
brid BnA configuration. J1331-I-305 was used as a 
flux calibrator, and J0958+655 was observed every 
15 minutes to calibrate the interferometer phases. 
Standard procedures for calibration, editing and 
imaging were carried out using the Astronomy Im- 
age Processing Software (AIPS) of NRAO. 

Images of M82 were produced for each obser- 
vation day as well as a composite image produced 
from the combination of all data, shown in Fig. 8. 
In order to study the compact sources in M82 at 
8.5 GHz, the extended emission larger than 4" was 



filtered out by limiting the corresponding (u,v) 
data range to baselines longer than 50 kA. The 
resulting resolutions of all images are subarcsec- 
ond (approximately 0.6" x 0.3"). 

We detected a low signal-to-noise radio source 
in two of the four observations at a position which 
is 1.2" from the position derived for X41.4+60 in 
our TOO observation, but only 0.7" from the po- 
sition quoted for X41.4+60 in Kaaret et al. (2001). 
Given the astrometric uncertainty of Chandra, we 
conclude that the source is a possible radio coun- 
terpart to X41. 4-1-60. This radio source is unre- 
solved in all observations and has a position of 
(J2000): 09 55 50.33, 69 40 46.3. The 8.5 GHz 
flux density of this source is 0.9 ± 0.1 mJy on 1 
Feb 2005 and 0.6 ± 0.1 mJy on 3 Feb 2005. The 
source is not detected on 29 Jan 2005 with an up- 
per limit of 0.3±0.1 mJy and a marginal detection 
of 0.5 ±0.1 mJy on 5 Feb 2005. Rodriguez- Rico et 
al. (2004) report the only other detection of a radio 
source at this position, 41.62+59.9 (J1950 coordi- 
nates name), with a flux density of 1.6 ±0.1 mJy at 
8.5 GHz. The source is 0.8" from the radio tran- 
sient reported by Kronberg & Wilkinson (1985). 

Radio emission is also detected from a source 
which falls within the 0.7" error radius of X42.3+59. 
This radio source, known as 42.21+59.0, has been 
detected previously by a number of radio studies of 
compact sources in M82. We measure the flux den- 
sity of 42.21+59.0 to be 2.9 pmO.l mJy during the 
flrst three observations and 2.4 + 0.1 mJy on 5 Feb 
2005. Values for the 8.5 GHz flux density of this 
source range from 1.9 mJy (Huang et al. 1994) to 
6.7 mJy (Rodriguez-Rico et al. 2004). The range 
of values is most likely due to the differing inter- 
ferometric arrays used for observation (which are 
sensitive to different spatial scales of radio emis- 
sion) and the different techniques used to account 
for the extended background emission surrounding 
the source. Based on its flat or slightly positive 
spectrum (between 5 and 15 GHz), extended size, 
the presence of OH and water mascrs and radio re- 
combination line emission, 42.21+59.0 is believed 
to be a giant HII region containing as many as 94 
05-type stars (Allen & Kronberg 1998; McDonald 
et al. 2002; Rodriguez- Rico et al. 2004). 



^The NRAO is a facility of the National Science Founda- 
tion operated under cooperative agreement by Associated 
Universities, Inc. 
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Fig. 8. — VLA 8.5 GHz image of a central region in M82 containing the two X-ray brightest sources (indicated 
by crosses of length 0.7", corresponding to Chandra positional errors). There is a possible (but low signal- 
to-noise) radio source at a displacement of 1.2" from the brightest X-ray source, X41. 4-1-60. The flux of 
this radio source is less than 1 mJy. The second brightest X-ray source, X42. 3-1-59, has a position which is 
coincident with the radio source 42.21-1-59.0 (its flux is ^^2.5 mJy and detected here with a signal-to-noise 
of 30). Contours show radio emission at levels of -3, 3, 5, 9, 12, 15, 20, 25, 35, 45, 60 and 80 times the rms 
level of 0.08 mJy/beam. The resolution of this image is 0.62" x 0.28", PA=79.8°. 



4. Discussion 
4.1. X-Ray Flare 

The fluxes measured by RXTE for the obser- 
vations before and after the Chandra observation 
are consistent with 3.2 x 10~^^ ergcm"^ s~^ in the 
2-10 keV band. The total flux in the entire Chan- 
dra image is 3.0 x lO^^^ergcm^^ s^^ in the 2- 
10 keV band, which is consistent within the er- 
rors on the flux estimate and the relative cal- 
ibration uncertainty between RXTE and Chan- 
dra. Of the total Chandra flux in this band, 
50% arises from X41.4+60 and 23% arises from 
X42.3+59. The peak flux measured by RXTE 
during the flare was 6.3 x 10^^^ ergcm^^ s^^ on 
MJD 53400.2. Assuming that a flux of 1.5 x 
10~^^ ergcni^^ s^^ arises from the remainder of 
the galaxy, the flux of the point source producing 
the flare would be 4.8 x lO^^^ergcm"^ s"^, equiv- 
alent to a isotropic luminosity of 7.6 x 10'*" ergs^* 
at a distance of 3.63 Mpc. This is consistent with 
the highest fluxes measured with Chandra from 
X41.4+60 (Kaaret et al. 2001). It appears hkely 
that X41.4+60 produced the flare observed with 
RXTE. 



The energy spectrum that we find for X41.4+60 
is best described by a single absorbed power-law. 
We find no evidence for thermal emission from a 
disk or for Fe-K line emission. The spectrum, a 
featureless power-law with a relatively hard pho- 
ton index, is reminiscent of the spectra of blazars, 
but the low level of radio emission, even during the 
X-ray fiare, rules out a blazar identification for the 
source. Adopting the maximum fiux that we ob- 
serve for the radio source 41.62-1-59.9 as an upper 
limit on the radio fiux of X41. 44-60, we find that 
the X-ray to radio flux ratio is more than an or- 
der of magnitude larger than found for any known 
blazar (Landt et al. 2001). The absence of strong 
radio emission during the flare argues against it 
being a relativistically beamed jet source (Kording 
et al. 2002). 

It has been suggested that ULXs may be me- 
chanically beamed (King et al. 2001) or radiating 
at super-Eddington rates (Begelman 2002). The 
fiux from X41.4+60 in the 2-10 keV band at the 
peak of from fiare corresponds to a luminosity of 
7.6 X 10'*" erg s~^ at a distance of 3.63 Mpc, equiv- 
alent to the Eddington limit for a 580 Mq object. 
We note that this is the observed (absorbed) fiux. 
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Using the best fit spectral model from the Chan- 
dra observation, this 2-10 keV flux corresponds 
to an unabsorbcd flux of 8.7 x 10""'^"'^ ergcm"^ s~^ 
in the 0.3-10 keV band for a luminosity of 1.4 x 
lO^^ergs"-^ or 1.3 x 10~-^°ergcm~^ s~^ in the 0.1- 
20 keV band for a luminosity of 2.1 x 10"*^ ergs~^. 
The latter luminosities corresponds to the Edding- 
ton limits for 1060 Mq and 1600 Mq objects, re- 
spectively. Given that the masses of known stellar- 
mass black holes are all below 20 M© (McClin- 
tock & Remillard 2003), this would require the 
inferred luminosity to exceed Eddington by a fac- 
tor of least 30 and more likely 50-80. However, 
models of 'funnel' shaped disks, as required for 
mechanical beaming, indicate that the maximum 
flux enhancement is only a factor of « 5 (Misra 
& Sriaram 2003). This falls far short of the en- 
hancement needed to explain X41. 4-1-60 as a stan- 
dard black hole. The formation of photon bub- 
bles in magnetized, radiation pressure dominated 
disks may also lead to super-Eddington luminosi- 
ties. Originally, the potential to reach luminosi- 
ties up to a factor of 100 in excess of Eddington 
were suggested (Begelman 2002). However, the 
best currently available calculations indicate en- 
hancement factors of a few (Ohsuga et al. 2005; 
Watarai et al. 2005) or less than 10 (Ruszkowski 
& Begelman 2003). 

The question of whether or not black holes 
can radiate at substantially super-Eddington lu- 
minosities can also be addressed by examination 
of Galactic X-Ray binaries and active galactic nu- 
clei (AGN). CoUin & Kawaguchi (2004) find that 
the bolomctric luminosities of AGN saturate at a 
few times the Eddington luminosity. In a sample 
of 578 type 1 AGN, Warner et al. (2004) find that 
only the most extreme objects have Eddington ra- 
tios approaching 10. We note that errors in black 
hole mass measurement and the AGN distances 
may contribute to the highest Eddington ratios 
quoted. For Galactic X-ray binaries with mea- 
sured black hole masses, the highest luminosities 
observed are only a few times Eddington (McGlin- 
tock & Remillard 2003). It seems that the maxi- 
mum luminosities produced by accretion disks in 
nature exceed the Eddington limit generally by no 
more than a factor a few and in the most extreme 
objects by less than a factor of 10. 

Given these considerations, it appears unlikely 
that X41.4+60 can be interpreted as a standard 



stellar-mass black hole. It's extremely high lumi- 
nosity and strong variability identify it as, per- 
haps, the best intermediate-mass black hole can- 
didate amongst the ULXs. 

4.2. X42.3+59 

Our Chandra observation found X42.3+59 at 
the highest flux yet seen. The 0.3-10 keV flux 
corresponds to an isotropic luminosity of 1.1 x 
10^" erg s~^, a dramatic increase in brightness over 
what was seen in previous observations. The 
total 2-10 keV flux found from the central re- 
gion of M82 during the XMM-Newton observation 
where Strohmayer & Mushotzky (2003) detected 
a QPO was 2.1 x 10^^^ ergcm~^s~^. It is pos- 
sible that X42.3+59 contributed one third of the 
total XMM-Newton flux and was the origin of the 
QPOs. The required modulation would be 30%, 
which is comparable to rms fractions measured 
for low frequency QPOs, in the range 0.6-2 Hz, 
from GRS 1915-^105 (Morgan et al. 1997). We 
note that our Chandra timing power spectrum of 
X42.3+59 shows significant variability near 1 mHz 
and a possible QPO at 70 mHz, similar to that 
seen with XMM-Newton, while the timing noise 
from X41. 4-1-60 is consistent with Poisson fluctu- 
ations in this frequency range. This may suggest 
that X42.3-h59 and not X41.4-h60 is the source 
of the QPOs detected with XMM-Newton and 
RXTE. 

If the radiation is isotropic, then the minimum 
black hole mass for X42.3+59 is 85 Mq. Black 
hole X-ray binaries tend to produce QPOs while 
in the low/hard state, which corresponds to lu- 
minosities below 10% of the Eddington limit. If 
X42. 3-1-59 is the source of the QPO detected from 
the central region of M82 and was is at a simi- 
lar luminosity while the QPO is being produced, 
then the estimated black hole mass is of order 
1000 Mq. This may indicate the presence of a 
second intermediate-mass black hole candidate in 
M82. 

4.3. X-Ray Periodicity 

We have found an apparent 62 day periodicity 
in the X-ray light curve of M82, which likely arises 
from X41. 4-1-60, the brightest X-ray source in the 
galaxy. This periodicity could be due to orbital 
modulation, superorbital modulation, or an aperi- 
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odic phenomena. 

The X-ray emission of several black hole X-ray 
binaries is modulated at the orbital period. In par- 
ticular, Cygnus X-3 presents a clear case, known 
since 1980 (Eisner et al. 1980), of a modulation in 
X-rays at the orbital period. More recently, Boyd, 
Smale, Dolan (2001) found a modulation at the 
1.7 day orbital period of LMC X-3, a black hole 
X-ray binary with a high mass companion. Also, 
recent results on systems with giant companion 
stars suggest that such modulations are more com- 
mon than in low-mass systems. Smith, Heindl, 
& Swank (2002) found modulations at 12 days 
for IE 1740.7-2942 for 18 days for GRS 1758-258 
which they interpret as orbital (they also find su- 
perorbital modulations at around 600 days), and 
Corbet (2003) found a 24.06 day period in X-rays 
from GX 13+1, which he interprets as an orbital 
modulation. The modulations from these systems 
are roughly sinusoidal, but tend to have sharper 
peaks and troughs than pure sinusoids. The range 
of shapes is consistent with the shape of the mod- 
ulation that we observe from M82. X-ray modu- 
lation at the orbital period appears common for 
systems with giant companions. 

There are so-called super orbital modulations in 
the X-ray emission from X-ray binaries. For neu- 
tron star systems, the super orbital periods include 
the range around 62 days. However, to produce 
the flux observed from X41.4-f60 from an accret- 
ing neutron star would required a beaming factor 
of larger than 300. This is much larger than has 
been suggested to be theoretically possible in me- 
chanically beamed models, see above, and would 
require a beam with an angular spread of less than 
0.2°. For black hole X-ray binaries, the shortest 
superorbital period is 162 days (from SS 433 for 
which there is no dynamical mass measurement 
confirming its black hole status) and the super- 
orbital periods are more typically in the 300-600 
day range (for Cyg X-1, IE 1740.7-2942, and GRS 
1758-258). Therefore, the observed 62 day period 
lies outside the range of observed superorbital pe- 
riods for black hole X-ray binaries. We cannot 
exclude the possibility that the process produc- 
ing the 62 day modulation could be aperiodic. 
However, we are not aware of any process which 
would produce aperiodic, but relatively coherent, 
behavior with the time scale and modulation am- 
plitude seen. We conclude that the periodicity in 



the RXTE light curve likely indicates the orbital 
period of X41.4-h60. 

Possible orbital periodicities have been re- 
ported for a few other ULXs. The previous 
best candidate is CG X-1 in the Circinus galaxy 
(Bauer et al. 2001) for which (almost) three pe- 
riods of 7.5 hours were detected in an observa- 
tion of 16.7 hours. The period was confirmed in 
additional (shorter) Chandra observations and a 
BeppoSAX observation of 14.4 hours (Bianchi et 
al. 2002), and also in an XMM-Newton observa- 
tion of 28.9 hours (Weisskopf et al. 2004). The 
period is consistent with a main sequence com- 
panion with a mass near 1 Mq or a similarly low 
mass core He-burning companion. The light curve 
of CG X-1 resembles that of an AM Herculis type 
cataclysmic variable star, with sharp transitions 
between high and low fiux phases, and the source 
lies relatively close to the plane of the Milky Way 
b = —3.8°, so it may be a foreground AM Her 
system Weisskopf et al. (2004). 

Other reported detections of orbital modula- 
tions from ULXs are of substantially lower signifi- 
cance. Sugiho et al. (2001) report possible 31 hr or 
41 hr periods from a 154 hr ASCA observation of 
IC 342 X-1, but the signal is relatively weak and 
even extraction of the best period is ambiguous. 
Liu et al. (2002) report a possible 2.1 hr period 
from M51 X-7, but only two cycles were detected 
in a single Chandra observation of 4.2 hr. Dewan- 
gan et al. (2005) report a periodicity of 1.6 hr from 
the same M51 source in a 5.8 hr XMM-Newton 
observation obtained 2.5 years after the Chandra 
observation, so it appears unlikely that the mod- 
ulation is related to the binary orbit and, instead, 
is either a quasiperiodic or red-noise modulation 
of the flux. David et al. (2005) report a suggestion 
of a period of 8-10 hr from a ULX in the ellipti- 
cal galaxy NGC 3379, but this is based on a 9 hr 
observation - a single orbital period. Quasiperi- 
ods have been seen from several ULXs including 
a ~ 2 hr quasiperiod from a source in NGC 7714 
(Soria & Motch 2004) and a ^ 1 — 2 hr quasiperiod 
from a source in M74 (Krauss et al. 2005). How- 
ever, these quasiperiods are not stable enough to 
be interpreted as orbital modulations. 

We note that very few ULXs have been moni- 
tored in a manner which would permit detection of 
modulations on time scales of tens of days. Indeed, 
the program described here is the only such ULX 
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monitoring of which we are aware. Therefore, the 
lack of detections of modulations on time scales of 
tens of days may simply be an observational ar- 
tifact. Monitoring programs extending over 100- 
300 days with observations every day or few days 
would be of great interest. We note that the 
best available optical constraint on the nature of 
a ULX companion star comes from an HST spec- 
trum from which Liu, Brcgman, & Seitzer (2004) 
identify the companion to NGC 5204 X-1 as a BO 
lb supergiant. They predict an orbital period on 
the order of 10 days. This provides additional mo- 
tivation for searching for X-ray modulations from 
ULXs with time scales of tens of days. 

4.4. Nature of the ULX companion star 
and formation of the ULX 

We note that X41.4+60 is hkely a Roche-lobe 
overflow system because wind-fed accretion can- 
not provide a sufficiently high mass accretion rate 
(Kaaret et al. 2004). In this case, the orbital 
period directly constrains the average density of 
the companion star according to the equation 
p ~ 115(P/hours)^^ gcm^"^ (Frank, King & Raine 
2002). If the orbital period is 62 days, then 
the mean density of the companion star is 5 x 
10~^gcm^'^. This excludes main sequence stars, 
but would be compatible with giant stars or super- 
giant stars. If, instead, the 125.7 day periodicity is 
the orbital period, then the companion star den- 
sity would be 1.3xl0^^gcm^^ and the companion 
star must, again, be a giant or supergiant star. 

As noted above, the only other ULX with a rea- 
sonably firm identification of the spectral type of 
the companion is NGC 5204 X-1 for which (Liu, 
Bregman, & Seitzer 2004) find a spectral type of 
BO lb based on a far-ultraviolet spectrum. This 
may suggest that giant or supergiant compan- 
ion stars are common in ULXs. However, more 
companion spectral types must be determined to 
drawn any firm conclusions. The spectral types 
for the companions to some other ULXs are been 
estimated using optical colors, but these classifi- 
cations are suspect because optical emission from 
the accretion disk could affect the colors (Kaaret 
2005). The density of 5 x 10~^gcm~^ would be 
consistent with standard late K giant or late A 
supergiant stars. However, the properties of the 
companion star are significantly affected by the 
high rate of mass transfer, and stellar evolution 



calculations taking into account the mass transfer 
are required to properly understand the evolution- 
ary state of the companion star. 

Stellar evolution calculations performed by Li 
(2004) for a star with an initial mass of 15 Mq 
orbiting a 1000 Mq black hole show that the mass 
transfer rate exceeds 10"'* M© yr~^ when the com- 
panion evolves through the giant phase at 10 Myr 
after formation of the binary. Similar results are 
obtained by Portcgics Zwart et al. (2004b). This 
mass transfer rate is more than sufficient to fuel 
the observed X-ray luminosity of 10** ergs~*. The 
lifetime of the giant phase is relatively short, on 
the order of 10^ years. This may suggest that 
the ULX phase is rather short lived and that sys- 
tems with very high luminosities, on the order 
of 10**ergs~*, are rare (Kaaret, Simct, & Lang 
2006). 

As noted in the introduction, X41.4+60 lies 
near and possibly within the super star cluster 
MGG 11. McCrady, Gilbert, & Graham (2003) 
estimate that MGG 11 is 7-12 Myr. Hence, the 
time required for the companion to evolve to the 
giant phase is consistent with the age of the clus- 
ter. 

A possible mechanism for the formation of in- 
termediate mass black holes with masses greater 
than IOOMq is stellar collisions in the cores of 
dense stellar clusters (Portegies Zwart et al. 1999; 
Taniguchi et al. 2000; Miller & Hamilton 2002). 
Portegies Zwart et al. (2004a) found that the ex- 
tremely compact size of MGG 11, the half light ra- 
dius is only 1.2 pc, causes massive stars to rapidly 
sink to the cluster center via dynamical friction. 
The time for sinking is shorter than the lifetimes 
of the stars, so a large number of massive stars ac- 
cumulate in the cluster center where runaway col- 
lisions produce a star of ^ 1OOOM0 that collapses 
to an intermediate mass black hole. This process 
does not occur in less dense clusters. The coin- 
cidence of the X41. 4-1-60 with a super star clus- 
ter with the correct properties for the production 
of an intermediate mass black holes favors run- 
away collisions as the production mechanism of the 
ULX. Other models for intermediate mass black 
hole formation, such as remnants of the first gen- 
eration of stars (Madau & Rees 2001) or as the nu- 
clei of captured satellite galaxies (King & Dehnen 
2004), would not predict such a coincidence. 
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